Dynamics beyond completely positive maps: some properties and applications 
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Maps that are not completely positive (CP) are often useful to describe the dynamics of open 
systems. An apparent violation of complete positivity can occur because there are prior correlations 
of the principal system with the environment, or if the applied transformation is correlated with the 
state of the system. We provide a physically motivated definition of accessible non-CP maps and 
derive two necessary conditions for a map to be accessible. We also show that entanglement between 
the system and the environment is not necessary to generate a non-CP dynamics. We describe two 
simple approximations that may be sufficient for some problems in process tomography, and then 
outline what these methods may be able to teli us in other situations where non-CP dynamics 
naturally arise. 



I. INTRODUCTION 

All real world systems interact to some extent with 
their environments, so they are said to be "open" [Ï],Q,Ü|. 
When the initial correlations with the environment can 
be neglected the evolution is well-described by a com- 
pletely positive map (CP-map). A CP-map can always 
be written in the Kraus form [J], 



P 



(D 



When the Kraus operators M a satisfy the completeness 
rclation, 



^MlM a = t, 



(2) 



the map is also trace preserving, so that tr p' = tr p. 

Howcvcr, if the system and the environment are ini- 
tially correlated the resulting reduced dynamics may not 
be CP SBSBISIS· Positive but non-CP maps also 
play an important role in characterizing the phenomenon 
of quantum entanglement [ni E2, EU, EU ■ Our goal is to 
investigate under what circumstances non-CP maps can 
describe an actual quantum dynamics and when (and if) 
the deviations from CP dynamics can be ignored. 

To simplify the exposition we consider finitc- 
dimensional systems. The combined state of a system 
(A) and its environment (B) can be represented in the 
Fano form [P5| 



tab 



— ( 



+ 



dAds 



AB 



<j*-A 



(3) 



Here the af, i = 1, d\ represent generators of SUÍcU), 
while the real vector ex of size d\ — 1 is the generalized 



Bloch vector of the reduced density operator p — tr bt~ab- 
Analogously, the af represent generators of SU(ds) and 
[3 of size d B — 1 denotes the generalized Bloch vector of the 
reduced density operator lu — ívatab- The correlations 
between subsystems A and B are characterized by [Ï6| 



(lij ~ ai8 3 )/dAdi 



(4) 



We assume that the overall evolution of tab is unitary. 
To specify a non-unitary dynamics of the system A we 
need to describe how it is embedded into a larger system 
AB. This is described by a map Ey such that, 



Ev(p) 



tr B (r) 



(5) 



which is called an assignment Q or an extension map 
A tensor product assignment with a an initial ui 
that is independent of p, on the auxiliary Hilbert space 
ri b followed by a unitary U leads to a CP-map [J] , 



p' = tr B {U(p®Lü°)U 1i ) 



(6) 



If the initial state of the environment lo is related to the 
initial state of the system p, then the reduced evolution 
of the system may be non-linear. For example [Ï7j . p 
may be an improper mixture (i.e., obtained by taking a 
partial trace from some larger system). This procedure 
yields an ensemble of pure states that depend on some 
classical parameter c, iPa(c) (g> ^b{c). Alternatively, the 
marginal state of the environment may be independent of 
p, but if the applied transformation $ depends on some 
parameter c, the final density matrix p' will not generally 
be equal to the result of applying the average of $ over 
c to p. This latter situation has arisen in the process 
tomography of a nuclear magnètic resonance quantum 
information processor (Ï8j . 

Part of the controversy surrounding non-CP maps in 
literature Q can be traced to ambiguities in the defini- 
tions of the extension maps. Moreover, the presence of 
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correlations may blur the boundary between the system 
and its environment. The main aim of this work is to 
introduce a class of non-CP maps that may be useful in 
the description of the dynamics of open systems corre- 
lated with the environment and to analyze some of their 
properties. 

The rest of this article is organized as follows. In the 
next section we provide a short review of some properties 
of positive linear maps. The notion of accessible maps is 
introduced in Section IIIII while their properties are in- 
vestigated in Section IIVI Some implications for process 
tomography are presented in Section [V] and a few other 
applications of non-CP maps are discussed in Section [VÏJ 

II. MAPS AND DYNAMICAL MATRICES 

In this section we recali several properties of lin- 
ear maps on the set of density matrices. A linear, 
hermiticity-preserving transformation <í> acting in the 
space of density matrices may be represented by the dy- 
namical (Choi) matrix -D($) |19l. 120] . 

Prnn = D ms . nt p st , (7) 

which has a number of useful properties [24| . The trace 
preserving condition (J2]) is equivalent to the following 
constraint on the partial trace of the dynamical matrix, 

D mS ;mt = Ò~st j (8) 

i n 

which implies A Q = Úa- Moreover, if the map is uni- 
tal, i.e., it maps the maximally mixed state into the max- 
imally mixed state, then 

s 

The dynamical matrix is Hermitian, {D^) ms . nt = 
D^t-ms: ancl due to a theorem of Choi [2(| its positivity is 
equivalent to the complete positivity of This property 
follows from the eigen-decomposition of the Choi matrix, 

D m s-nt = ^\a{M a ) ms (MÍ) tn , (10) 
a 

in which all the eigenvalues A a are non-negative. 

If the dynamical matrix D is not positive we can split 
its spectrum into positive and negative parts. This step 
allows us to represent a linear non-CP map as the differ- 
ence of two CP maps, called the difference form [9l [20|. 

p' = A + (p) - A_(p) (11) 
= x * M «P M l + J2 KM a pMl, (12) 

A a >0 A a <0 

where the maps A± are completely positive. 



Completely positive maps have another important 
property: distinguishability of the set of states 
(pi,p2, ■ ■ ■) does not improve under any CP map [23[. 
For example, the trace distance between two density ma- 
trices does not increase under any CP map A , \px — p% | 1 = 
tr|pi - p 2 | <tr|A(pi - P2 )\. 

Despite its definition as a mathematical tool, some 
matrix elements of D have a direct experimental signif- 
icance. For example, by a narrow-band làser resonant 
transition Z <-* Z' in a nitrogen-vacancy (NV) center in 
diamond, the fluorescence intensity autocorrelation func- 
tion 5 (2) (í) is 3 (2) (í) = D zz ,. zz ,/p z , z ,, where p x > z > is the 
steady-state population of the fluorescent substate Z' 

III. MAPS AND EMBEDDINGS 

A possibly non-linear and non-positive map $ that de- 
scribes a state transformation on A may be considered 
physically accessible if $(p) = tr bU ~tU\ with the embed- 
ding described by some E\>. The domain V of $ should 
be a finite volume (i.e., non-empty) subset of the set of 
all states of A, which we will call AÍa- This is the first 
non-trivial requirement: physically relevant maps should 
be identifiable by process tomography, and convex combi- 
nations of a tomographically complete set of states must 
span a finite- volume region of Ma- In addition, a non- 
positive map $ may only be accessible for states in its 
domain of positivity, where $(p) > 0. 

However, this definition is still too broad; some more 
restrictions on the assignment maps are essential. With- 
out these further conditions the definition of accessibil- 
ity remains trivial: any map becomes accessible on its 
domain of positivity. For example, the positive non- 
CP transposition map T : p 1— > p T , results from the 
extension E(p) — p (g) p T followed by the SWAP gate, 

Í^SWAp(p ® w )^swap = 10 ® P- ^ n arbitrary non-linear 
map p 1— > p final C an be realized by E{p) — p ® p final and 
the unitary SWAP. 

It might be objected that this construction is rather 
contrived. However, when we are setting the dials of 
our preparation apparatus in order to produce our to- 
mographically complete set of input states, we are giving 
that apparatus a complete, classical description of the 
state we would like it to produce. Once we have done 
that, there is no a priori reason why the apparatus should 
not produce extra copies of the state, or ones that have 
undergone some peculiar map. 

There is another way of demonstrating this point that 
draws on quantum information about the tomographi- 
cally complete set of input states instead of the classical 
description of the state invoked above. Suppose instead 
that we are given that the environment consists of in- 
finitely many copies of the state of the system (in a ten- 
sor product). Then the environment contains a complete 
(i.e., classical) description of the state, and can be used 
to implement an arbitrary map. This can be seen by not- 
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ing that we could use the copies in the environment to do 
exact state tomography and then proceed as above [22| . 
There is also a direct equivalence between partial quan- 
tum information about a state (in the form of finitely 
many duplicate copies of that state) and incomplete or 
"fuzzy" classical descriptions of that state, via optimal 
state tomography. 

Therefore we can implement any map that depends on 
detailed (classical) knowledge of the state provided we 
have access to an environment that contains a sufficiently 
large number of copies of that state. If the map requires 
knowledge of the state that is infmitely precise (i.e., in 
order to perform this map correctly, we must pick out 
a lower-dimensional subset of the set of density matri- 
ces) then an infmitc number of copies in the environment 
would be required. For example, if we know what p is ex- 
actly, we can always perform the map p i— > p T by another 
route: Find the unitary U such that UpW is diagonal, 
and then do UUpWU* = p T . However, this map depends 
on p, via U, and it requires exact knowledge of the eigen- 
states of p. There are also other ways of using múltiple 
copies of a state to perform exòtic maps on that state. 

We can draw a number of conclusions from these 
observations: 

(i) if the ancilla system lo depends on p, a non-CP map 
may arise from a tensor product assignment; 

(ii) the domain V may be the entire set of states Ma] 

(iii) Given enough copies of p in the environment, any 
map can be performed [22j |. 

Clearly realistic systems will not have environments 
that contain so much information about the state that 
point (iii) becomes a completely unmanageable problem, 
but how much information is it reasonable to assume the 
environment might have about the state? The informa- 
tion known to the environment is part of the assignment 
map, thus our first task is to demarcate the set of physi- 
cally reasonable assignment maps. We will begin by con- 
sidering the simplest, "linear" case where the assignment 
map is linear; i.e., the marginal state of the environment 
uj is independent of p and the correlations between sys- 
tem and environment can only be seen in the density 
matrix of the combined system tab- We will then con- 
sider the slightly more general "affine" case, before dis- 
cussing what approximations might be possible for the 
"non-linear" cases where u> is allowed to depend on p 
directly. 



A. The linear case 

The simplest scenario is an initial value problem in 
which the time evolution of different initial states of the 
system is analyzed given that the initial state of the en- 
vironment and the system-environment correlations are 
independent of p. If the initial Bloch vector of the system 
is a, the (3j could still depend on the on. If this is not the 



case (i.e., /3j = bj, a constant) then we can write 

Pj = bj, "fij = gij + aibj. (13) 

where gij — dAds^ij- In general, Tij may depend on p 
as well as the c^, but for the case when <?jj is a constant 
matrix, 7 depends only on the ai, via the second term in 

m- 

Under the action of a unitary U on the extended 
system, a useful form of the reduced dynamics of p 
is obtained using the following procedure proposed by 
Stelmachovic and Buzek [7|. 

Decompose the extended density matrix tab uúo a 
simple tensor product and the remaining term, 

TAB = P® v + {tab - P® uj). (14) 

The direct product term yields a CP-map, as u> is inde- 
pendent of p so this term is in Stinespring form. Using 
the spectral decomposition u> — ^2 v Pv\v)(v\ an d writing 
the partial trace as 

ix B (Up®uUÍ) = J2(v\VpÏU\v} PA (v\VpÏUÏ\p), (15) 

where the \p) form an orthonormal basis in Tíb, and 
defming M M „ = {p\^/p~U\v), one obtains Eq. (JTJ) after 
merging the double index pv into a single index, a. Hence 
the affine form of $>(p) — p' is given [7] by 

$(p) = tr B [U(E v (p))rf] 

= £ M aP Ml + £ £ -rÍr^\U*f ® afU^). 



B. The affine case 

The assumption that the initial environmental 
marginal uj° cannot depend on the initial state p° is 
rather strong; for realistic (and poorly characterized) sys- 
tems, we should not rule out any functional dependenec 
a priori, thus we should treat lo — u>°(p ) unless we have 
good reason to assume the system does, in fact, behave 
like the linear case. 

As noted above, if the environment has access to in- 
fmitely many copies of p , then any dynamical map can 
be induced on the system and the problem is intractable. 
However, that scenario is also rather unnatural. The 
question now becomes: can we make some physically 
well-motivated assumption about this system that will 
also make the problem tractable? We will make the as- 
sumption that although the environment "knows" some- 
thing about the state, it only knows a little information. 
We believe that this is a plausible assumption to make for 
physically reasonable systems and we will proceed with 
our analysis on that basis. 

In practice the state tab may result from the evolution 
of p° £g) u>° under some (imperfectly controlled) unitary V 
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which acts on the combined system. This can be a more 
rcalistic description of scenarios such as a gate being ap- 
plied to a qubit that was stored in an imperfect quantum 
memory than a simple CP map. When the desired gate is 
eventually applied, the target state p is only an approx- 
imation to the intended state p°. The dynamical matrix 
for the process p° <— > p is readily obtained as 



example above: p 



„ final 



Of course, those examples 



D, 



ac;bd — ^aa-c-fV^^gUJ^g. 



(17) 



Likewise, the evolution of the environment is described by 
a similar expression, with the environment 's dynamical 
matrix depending on p° : 



q,7,{5 



(18) 



In genèric cases D describes a one-to-one mapping 
p° i— ► p and therefore defines a unique transformation 
between the various coefficients. The vàlues of the coef- 
ficients for the transformed matrix p can be obtained by 
projection onto the original basis in the usual way: 

^ = trafp, fa = trafui, 7^ = tr(af <g> crfr), 

(19) 

Inverting the transformation in Eq. (|17p and using the 
dynamical matrix for the environment introduced in 
Eq. (fï8|) above, we obtain the effective assignment map 
p 1 ► tab with 



= y^-BjfcQ!fc, = E Gí -J feQ!fe ' (20) 

k 



with some constant coefficients Bj k and Gijk, which are 
independent of the a 2 : but do depend on /3q. 

If D is not a one-to-one map, then p does not determino 
lo uniquely. Since all the relationships between the state 
parameters are linear, the effective assignment map in 
this case is given by 



k 



7 Un 



■^GyfcÜt, (21) 



which is actually the most general form of an affine as- 
signment map, where the constants are subject to the 
positivity constraints, as before. 



C. More general, non-linear cases. . . 

An assignment map may also describe a preparation or 
be part of the specification of an imperfect quantum gate. 
In this case the correlation with the environment and/or 
its dependence on the state of the system are established 
either by design or accident, and the relations between p, 
lo and r are potentially arbitrary. Allowing an arbitrary 
assignment map (even of the second degree) results in 
declaring all maps physically accessible, in principle, in- 
cluding perfect cloning, in a similar way to the transpose 



do not exclude any functional forms of non-linear assign- 
ment maps. However, in practice the criteria for "rea- 
sonable" restrictions are rather ad hoc: while it is quite 
clear that the completely unconstrained case that allows 
perfect cloning should be excluded, the feasibility of less 
outlandish assignment maps is determined by the details 
of the actual physical system and the preparation meth- 
ods. We leave those cases (like NMR [l8| or NV-centre 
qubits [3Ç]]) for a later study. We discuss a particular ly 
important example in in Section IVl 

Our discussion of the assignment maps can be summa- 
rized by the following definition. 

Definition A map $ defined by an assignment Ey(p) 
and a unitary U, 

$(p)=tr B [U(E v (p))tf}, (22) 

is called affinely accessible if the assignment map satisfi.es 
the linearity conditions in Eq. (|2ip and the unitary U 
does not depend on the initial state p. 



IV. SOME PROPERTIES OF ACCESSIBLE 
MAPS 

The affine form of a general linearly accessible map 
from Eq. (|16[) can be written 

<D(p) = E M aP Ml + E E r ^ ® ^f^\p). 

(23) 

Since tr Ugí ® o-jW = tr <Tí £3 Cj — 0, the inhomogeneous 
part is traceless and can be expressed as £ • a A . The CP- 
map term is trace-preserving. Moreover, since the final 
state p' is Hermitian, all the coefficients are real. As a 
result, the dynamical matrix is given by 



(24) 



/ í. . V 



and the vector <f S R d ^ _1 can be obtained by a compar- 
ison with Eq. (fï6|) . 

If the building blocks of a non-CP map (namely, an 
cxtcnsion map and a unitary evolution of the combined 
system) are linear and independent of the state p, then 
the map $(p) is linear and the resulting dynamical matrix 
is independent of p. Indeed, it is easy to verify that 
$(cpi + (1 - c)p 2 ) = c$(pi) + (1 - c)$(p 2 ). However, 
the affine form may still depend on the initial state of 
the combined system via the correlation tensor T and/or 
lo and it will therefore appear non-linear (compare with 
[27]]). In particular, a general assignment map of our 
example above, in Eq. ([26]) would lead to a quadratic 
dependence of Ty on p, 



= — — — — (gij + '^ / {Gi ] k - Sikbj)ak - E Bjkonak). 

" V ' B k k 

(25) 
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Another example of an assignment map is the exten- 
sion tab of a onc-qubit mixed state, p = i(l + a • a) 
to 



1 

TAB = 1 



AB 



3 



The operator tab is positive if 

< \a\ < <w = U/4-3|a| 2 -l)/3. 



(26) 



(27) 



Thus for any fixed a, the domain V of the map £7 is equal 
to the ball of radius \a\ = y(ï + a)(l — 3a), centred on 
the maximally mixed state. 

If the matrix Í7 commutes with the correlation term, 
which in our example means [U, J2 t of ® af\ = 0, the 
resulting map is trivially completely positive: the inho- 
mogeneous part £ • a is zero (25|. 

Now consider the two-qubit unitary rotation, 



U 



'1 0\ 

cosé» sinó» 

-sia 5 cosfl 

K 1/ 



(28) 



The inhomogeneous part of the resulting <E> is now non- 
zero: 



£(M)-<7 



1 /asin2é> 

2 l 





-asin29 



(29) 



however that alone is not enough to show that the map is 
non-CP. A typical situation is presented in Fig. [TJ which 
shows the spcctrum of the dynamical matrix as a function 
of the phase 9. The map $ is not CP for some vàlues of 9; 
for 9 = 7r/4 the affine form of $ is merely an inconvenient 
way to write a CP map, whereas for 9 = tt the map is an 
accessible, genuinely non-CP transformation. 

The conventional wisdom is that non-CP maps can 
only happen if the system is initially entangled with the 
ancilla. Our example consists of a two-qubit system and 
it is possible to detect entan glem ent using the positive 
partial transposition criterion |ll|, |l2| . The results of this 
test show that the state r is always classically correlated 
for a > (For a < the state r may be entangled, e.g., 
the case a = 0, a=— 1/4 corresponds to the triplet state 
\f r+ ). Moreover, even if U(9) leads to a non-CP map, the 
state U(9)tII' í (9) is still unentangled. 

Proposition 1. The existence of the affine form as in 
Eq. (|16[) with a trace-preserving completely positive map 
(that has at most linear dependence on p) and a trace- 
less inhomogeneous part (that is at most quadratic, as in 
Eq. (I25|) ) is a necessary condition for $ to be accessible 
through a linear extension. 

For a finite-dimensional ancilla the result follows from 
the definition of accessibility and Eq. (f2"3"l) . In the infmite- 
dimensional case we decompose a density operator tab 
as in Eq. (|14p . After applying a unitary operation to tab 
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FIG. 1: Three of the four eigenvalues of the dynamical matrix 
D(8; a — 0.2) become negative. Each line is one eigenvalue 
for a fixed value of a = 0.2 and 6 £ [0, 2n]. 

and taking the partial trace, the first term on the RHS 
yields a trace-preserving CP-map, and the second one is 
traceless and inhomogeneous. □ 
The conditions for Proposition 1 are non-trivial, as 
the next example shows. Consider the transposition of a 
qubit. Since it is a non-CP map, it has a difference form, 

T(p) = (o- + pa + + a-pa- + a x pa x /2) - (io y )p(io y ) i < /2, 

(30) 

where o~± = (1 ± a z )/2. However, this map has no affine 
form and it is therefore not accessible, as we now demon- 
strate. 

According to Proposition 1 the dynamical matrix of 
Dt should be decomposable into a trace-preserving CP 
part L(Ç) = ^2 a M a (g> Ml and the inhomogeneous part 
£ • a ® ï, as in Eq. 



(31) 



All the eigenvalues of L(£) must be non-negative for 
T(p) to be linearly accessible. However, a direct cal- 
culation shows that there is a negative eigenvalue A_ = 
— (1 + Ç 2 ) 1 , so L is not positive for any £, and thus 
is not linearly accessible. A non-linear realization with 
E{p) = p <g> p T and U = SWAP can be written in the 
form p T = Yla=i M a (p) pM a (p)^ , but finding the matri- 
ces M a (p) requires a complete (i.e. classical) knowledge 
of the eigenbasis of p. 

Moreover, it is interesting to note that a map 



T 1 = pA dcp + (1 - p)T, 0<p<í, 



(32) 



where Ad cp (p) = 1/2 is the totally depolarizing channel, 
is linearly accessible only when it is actually CP, i.e., for 
2/3 <p < 1. 

Proposition 2. Any unital map that has a state- 
independent affine form is completely positive. 

Any such map <E> has a form of Eq. (fTü]) . Following [26[ 
we decompose A as 



A( /9 )=A (p) + çtA)·a ! ) 



(33) 
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where Ao is a unital part of the map and ç represents a 
translation of the generalized Bloch vector that depends 
only on the map A. The requirement <&(!) = 1 implics 
(ç + 1) • <t = 0, so $(/?) = Ao(p) and $ is a CP map. □ 
We note that this simple result can be immediately ap- 
plied in the quantum causal histories approach to quan- 
tum gravity to reduce the number of independent axioms 
that characterize the dynamics of subsystems [28[ . 

V. SOME IMPLICATIONS FOR PROCESS 
TOMOGRAPHY 

Imperfections in the preparation procedure may lead 
to non- linear correlations with the environment. If those 
correlations are sufficiently weak the assignment map 
may be written 

0j = /3/ n + e/3j(a), 7tí = 7l f +£7Í(a), £ « 1, (34) 

where f3j and 7^ are arbitrary functions of a. 

Consider now a weakly correlated and weakly interact- 
ing subsystems A and B. Then the embedding p 1— * r is 
given by Eq. (34]) with /3j in = /3° and 7/j 11 = a t f3°. As a 
result, the state of the environment is given by 

u) =u° + eui 1 (p), trw^O, (35) 

and the system-environment correlations are described 
(to first order) by 

T {j = eT\ 3 = -JL- (71. - ai p) ) . (36) 

When the system-environment interaction is weak, 

H AB =H A ® 1+riHint, 77 <1, (37) 

where to simplify the notation we ignore the self- 
Hamiltonian of the environment. We also assume that 
the Hamiltonians are state-independent. The most gen- 
eral form of the Hamiltonian is 

h a = y. H ^ = H h a a ? ® °f • ( 38 ) 

i ij 

It is worth noting that the assignment map p 1— ► r is 
non-linear at the first order of e. 

Proposition 3. At the first order of the perturbation 
theory in both the correlation strength and the interaction 
strength a reduced dynamics that results from the above 
assignment is linear and CP. 

The unitary time-evolution operator is 

U AB = exp(~it[H A ®1 + r)H int ]). (39) 

To leading order in r\ it becomes 

Uab = U a ® 1b(1ab + vQi + n 2 Q 2 + ...) 

= U A ®tB+vOi+if0 2 + ..., (40) 



with 

00 

ij fe=l 

where the coefficients Cy can be found using the Baker- 
Campbell-Hausdorff formula [3l| and SU(ds) commuta- 
tion relations. 

Assume for simplicity that p° B has maximal rank and 
is non-degenerate. The Kraus matrices in the affine form 
of the reduced dynamics of Eq. (fTïï|) are independent of 
Pa up to the second-order terms. Indeed, 

Mftv w (mKVpS + epl/2y/fi)(U A ® 1 + r,0i)\v), (42) 

where \p) and are the unperturbed eigenvectors and 
eigenvalues, respectively, of u>° and the indices pv label 
the set of matrices, not the entries of a matrix. To first 
order we have 

*V - yM'UA5 ti u + ^=U A 5 l·lv +r l ^MOM. (43) 

2 Vp" 

When we substitute this expression into Eq. (JTJ) , the 
terms that are first order in e will cancel; thus (to first 
order) we can write 

= M° u + tiy/fi(ii\Oi\v), (44) 

where the matrices M° form the Kraus representation 
of Ua- Next, the inhomogeneous part (as defined in 
Eq. (ÈSD) 

Í- à A =sJ2J2 T lM U ^f ® °fu\ B \v), (45) 

is zero at the first order, because 

Y J U A o?U\{p\of\p)=Q. (46) 

Hence at the first order of the perturbation theory the 
reduced evolution is still linear and CP. □ 

If A is a cluster of qubits, B is its environment and the 
reduced dynamics represents a physical realization of the 
perfect gate Ua, a high fidelity (of the actual outputs p° A nt 
with respect with to the ideal outputs Ua^a)) a U° ws us 
to conclude that the first order perturbative expansion is 
vàlid. Hence Proposition 3 applies and the gate should 
be described by a CP map. 

The raw tomographic data often yield non-positive dy- 
namical matrices, which are usually considered unphysi- 
cal [29|, [3(| ■ A maximum-likelihood estimation or other 
such technique is used to convert the experimental data 
into a (positive) dynamical matrix [22|, |29j . We see that 
this can be justified for characterizing actual high-fidelity 
implementations of "known" gates. However, when e and 
77 cannot be considered "small" a different template (such 
as a difference form) should be used to fit the data when 
attempting linear inversion process tomography. 
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VI. OTHER APPLICATIONS 

In this section we will examine a couple of other appli- 
cations of induced dynamical maps that are non-CP. 

A. Dynamical decoupling 

The preservation of quantum memory by dynamical 
decoupling from the environment [32I ] clearly indicates 
that the reduced dynamics is non-CP. Let us revisit a 
simplified description of dynamical decoupling. Consider 
a system and the environment, initially in the state p®u>. 
Let the interaction Hamiltonian be Hab 1 so the evolution 
in the interaction picture is given by U = exp(— íHab^)- 
Assume that it is possible to produce a unitary map (such 
as an NMR pulse) P = Pa <8> 1b that anti-commutes with 
Hab, 

PH AB P ] = -H AB - (47) 

For example, if Hab = guf <8> Ob, where Ob is some 
operator that only acts on the environment and g is a 
coupling constant, then Pa = cr x ■ The pulse sequence 

p e -iH AB t p ï = jji ( 4g ) 

will reverse the original evolution U provided that the 
pulse P has negligible duration. Hence we will obtain 
p(2t) = p(t = 0). 

From the point of view of the system A the above evo- 
lution appears to be an accessible (and possibly state- 
dependent) non-CP map. Let us assume for simplic- 
ity that the system is a single qubit parametrized by 
the Bloch vector, the environment is finite-dimcnsional 
and was originally in a completely mixed state. We will 
also assume that from the point of view of the envi- 
ronment alone, the evolution U is a unital CP map, so 
u>(t) = lü = t/ds- Then the evolution under U of the 
initial state of p® l/ds leads to 

TAB{t) = ^{tAB+Y. aíS f^t® tB + Y. a ^% a t® a k), 

(49) 

where 

U(a£ ^)C/t = <>n ® <■ (50) 

It is assumed here that the coefhcients s are time- 
dependent, the indices can take arbitrary integer vàlues, 
p = 0,1,2, ... , and j = 1, 2, . . . , while the identity is de- 
noted by do : = 1- It is easy to see that the time-reversed 
evolution is given by a generalized inverse, 

U i {<rf®<r*)U = c%a£®<rï, where = 5% 6$. 

.( 51 ) 

The state tab (t) can be thought of as an image of a linear 
assignment map that was applied to 

p(t) = ^(l A +a i s%af) (52) 



The pulses of Eq. (|48]) result in a non-CP map $ that 
restores the state p(t = 0). 

B. Quantum channels 

A noisy quantum channel is usually modelled as a 
trace-preserving completely positive map. The informa- 
tion transmission from A(lice) to B(ob) can be then rep- 
resented as an isometry between Alice's Hilbert space and 
the Hilbert spaces of Bob and the environment, 

V :Ha^Hb®H c , (53) 

that is followed by a partial trace over the Hilbert space 
TLc (which is controlled by their colleague, Charlie). Sev- 
er al different channel capacities are defined, depending 
on the types of information and the resources that are 
at the disposal of the communicating parties. A typical 
message that is represented by a pure state ip € TL is 
block-encoded by Alice (with a block size n) through the 
operation A : 7í — > and is then sent through the 

channel, V® n . 

Recently there has been some interest in the capacities 
of channels that are assisted by a "friendly" environment, 
which can measure states on Tic and communicate the re- 
sult to Bob, thus potentially increasing one of the channel 
capacities [33]. In these scenarios Charlie (who observes 
the environment) measures the environment before Bob 
attempts to recover the information. The measurement 
is described by a POVM E x , 

E x >0, J2 E * = 1 c n > ( 54 ) 

X 

on 7í®" and the outeome x is communicated to Bob. The 
latter acts with the map 1Z X on his output, so the overall 
state transformation is given by 

*M = ^K*(tra»\y 9n A(1>)VÍ* n (lf l ®E x )]) . (55) 

X 

Such improvements in the distinguishability (and 
hence the capacity) show that from the point of view 
of the reduced states on Hb, the overall procedure that 
starts with Charlie's measurement and ends with 3>(ip) 
must be non-CP. 

In certain situations the process of encoding the "ideal" 
states (e.g., qubits) by Alice into the physical carriers 
may involve additional degrees of freedom. For example, 
when qubits are realized as a photon's polarization and 
the finite size and spread of the wave packets is taken into 
account, the operations A and V cannot be separated [35| 
and the channel is not described by a CP map. Even in 
the absence of other sources of noise, continuous degrees 
of freedom can play the role of an environment, while a 
subsequent passage through lcnses may lcad to a non-CP 
evolution of the polarization degrees of freedom. 

VII. OPEN QÜESTIONS 

The structure and applications of non-CP maps mèrit 
further investigation, particularly for the analysis of re- 
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alistic quantum gates. While it was recently shown [36J 
that for certain classes of extensions to separable states 
the reduced dynamics is always CP, there are still sev- 
eral important open qüestions. What conditions are suf- 
ficient for a map to be linearly accessible? What is the 
structure of the set of all linearly accessible maps? Wc 
have seen that the improvements in the distinguishabil- 
ity of quantum signals when the parties communicate, 
improvements in channel capacity or the preservation of 
quantum memory by dynamical decoupling from the en- 
vironment are examples of non-CP maps. Their proper- 
ties should be investigated in detail. 

Another group of qüestions is related to following the 
reduced dynamics through time. A CP evolution forms 
a quantum dynamical semi-group and corresponds to a 
Lindblad-type master equation |H, @, It is still an 
open question how a non-Markovian evolution is linkcd 
to non-CP maps @,[37j]. 
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